READ IMAGE matlab code :
mpo663_hw3_read_satellite_image.png

Figure 2: FFT of figure 1
mpo663_hw3_fft_image_figure2.png

Figure 3: wavenumber versus power [abs(Ihatsh)] : I may want to change y-axis (unit) later.
mpo663_power_wavenumber_figure3.png
In the high wave-number region: slope seems to proportional to -3 (higher than 10^3 wavenumber)
Between (approximately) 10^2 and 10^3 : the slope is ~ wavenumber^-2
In the lower wave-number region (lower than 10^-2) : slope seems to proportional to -1.
??? where is -5/3 slope??? It is synoptic picture....so 2D?

%---------------------------------------------------------------------------------------------------------------
I corrected the above figure. But...still not sure (it looks a bit strange to me)
x-axis: (log)wavenumber
y-axis: wavenumber * (sum(Ihatsh^2))
In general, the data follows -5/3 slope between wavenumber 2-6 (Does it make sense? it's synoptic scale...).
mpo663_hw3_sil.png
------------------------------------------This is not for this hw-------------------------------------------------
Figure 2: General feature of wind
black dashed line : 0 longitude
red dashed line : 160W (index=81), longitude=200
mpo663_hw3_wind_general.png

Figure 5 (1D FFT) : Frequency versus power plot at 160W (longitude=200)
To get the time at each peak I used Time = 1/f
mpo663_hw3_1d_fft_fig5.png

Figure 3 : Time series of wind at 160W
mpo663_hw3_wind_at_160W.png